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Abstract 

We analyze stability conditions of "Maclaurin flows" (self-gravitating, barotropic, two di- 
mensional, stationary streams moving in closed loops around a point) by minimizing their 
energy, subject to fixing all the constants of the motion including mass and circulations. 
Necessary and sufficient conditions of stability are obtained when gyroscopic terms in the 
perturbed Lagrangian are zero. To illustrate and check the properties of this new energy 
principle, we have calculated the stability limits of an ordinary Maclaurin disk whose dy- 
namical stability limits are known. Perturbations are in the plane of the disk. We find 
all necessary and sufficient conditions of stability for single mode symmetrical or antisym- 
metrical perturbations. The limits of stability are identical with those given by a dynam- 
ical analysis. Regarding mixed types of perturbations the maximally constrained energy 
principle give for some the necessary and sufficient condition of stability, for others only 
sufficient conditions of stability. The application of the new energy principle to Maclaurin 
disks shows the method to be as powerful as the method of dynamical perturbations. 

Key words: Energy variational principle; Self-gravitating systems; Stability of fluids. 
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1. Introduction 

The main purpose of this work is to illustrate on a test case our principle of "maxi- 
mally constrained" energy minimum (Katz, Inagaki and Yahalom 1993), and find stability 
limits of stationnary motions of barotropic flows. Paper I dealt with three dimensional 
flows and gave no application. Here we consider two dimensional "Maclaurin Flows" and 
apply the method to Maclaurin disks. With "Maclaurin Flows" we mean some sort of 
generalized rotating streams, with closed loops around a point, but no special symmetry. 
The simplest and most symmetric example of a Maclaurin flow is a Maclaurin disk. Other 
two dimensional fluid models for galaxies are Maclaurin flows as well. Stationary flows 
are "stationary points" of their energy. Flows are stable if their energy is minimum, com- 
pared to that of all slightly different stationary trial conflgurations . By constraining trial 
conflgurations to satisfy all the constraints that a dynamical flow would satisfy: linear mo- 
mentum, angular momentum, mass conservation and circulation, we may greatly enhance 
the value of the energy principle. Energy principles give only sufficient conditions of sta- 
bility. However, we know that if the "gyroscopic term" in the Lagrangian of a dynamically 
perturbed configuration is zero, a minimum of energy becomes a necessary and sufficient 
condition f 

Our test model is a Maclaurin disk. This is one of those rare models (Binney & 
Tremaine 1987) were a detailed dynamical stability analysis is possible. It is thus the one 
valuable test case in which the details of our method and its results can be compared with 
a different, independant calculation. Both method use spherical harmonic decompositions 
of perturbations or {l,m) - modes. 
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What we find are the same stability limits for all symmetrical and antisymmetrical 
single modes as those given by the dynamical analysis of Binney and Tremaine. Other 
models have been tested which give also necessary and suflBcient conditions of stability. 
Thus, Yahalom (1993) found the stability limits of a uniform rotating sheet and of two 
dimensional Raleigh flows, by the present technique, to be the same as that of a dynam- 
ical analysis. The new energy principle appears thus as a powerful instrument that may 
complete and occasionally replace the method of linear perturbation analysis. This is the 
main point of this work. 

2. Description of Stationary Maclaurin Flows, and of the Energy Principle for 
their Stability. 

2.1. Lagrange Variables and Equations of Motion for Stationary Flows. 

Equations of motion in three and two dimensions are similar, though in two dimensions 
angular velocity, angular- momentum and vorticity have only one component, in the z 
direction, and there are less variables. We use the following notations: for the positions of 
fluid elements R or (x^) = {x, y); K, L, ... = 1, 2, the density of matter per unit surface is 
cr; W is the velocity field in inertial coordinates and U the velocity relative to coordinates 
moving with uniform velocity b and rotating with constant angular velocity Qc these are 
the coordinates in which configurations appear stationary. Thus, by definition. 



W = U + fi, 



(2.1.a) 



ffc = b + lie X -R 



(2.1.6) 



rig — l2;fic* 



(2.1.C) 
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The mass is conserved in moving coordinates 

V • (aU) = (2.2) 
On the boundary of free self-gravitating fluids the density and the pressure are zero 

a\B = 0, (2.3.a) 

P\b = (2.3.6) 

and in flows of astrophysical interest the velocity of sound is there zero also, 

dP 

a^Uo = 0. (2.4) 

Notice that (2.2) and (2.3) imply 

{U-Va)\B = 0. (2.5) 

Since V(j|b is normal to the boundary as follows from (2.3a), (2.5) means, as one expects, 
that U is tangent to the boundary. Because of (2. 3. a) an integral of a divergence of the 
form V • {(tF) must thus also be zero by virtue of Green's identity: 



V • {aF)d^x = f aF-dB = (2.6) 
s Jb 

provided F is single valued and continuous in S and on B. Euler's equations for stationary 
motions in moving coordinates are well known and can be written: 

U -VW + ncXW + Vih + ^) = (2.7) 

where the speciflc enthalpy 

h= [^ = e{a)-^, (2.8) 
J a a 
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£ is the specific internal energy. The gravitational potential $ is, in the plane of the motion, 
given by 

$(^) = -g\ -^^f^- (2-9) 

The vorticity of the flow 

cl; = rot = rot C/ + (2.10.a) 
a; = uA^ (2.10.6) 

2.2. Lagrange Variables for Maclaurin Flows. 

The labels of fluid elements are appropiate Lagrange variables in fluid mechanics 
(Lamb 1945). A continuous labelling of fluid elements depends on the particular topology 
of the flow. Here we describe a special labelling, and a representation of the velocity fleld 
for which mass conservation and circulation conservation hold automatically. We consider 
Maclaurin flows that move in closed loops around a flxed point. Suppose that, in some 

— * 

coordinates moving with constant velocity b and rotating with constant angular velocity 
Qc, we are given the density a{x^) and the relative velocity of the flow U{x^). We do not 
suppose that a and U represents a solution of Euler's equations (2.7). We rather regard cr 
and [/ as a trial conflguration of the flow. With C/, b and VLc we deflne a velocity W{x^) 
in inertial coordinates by equation (2.1). Given a and W we can construct a quantity 

A=- (2.11) 

where uj is the vorticity deflned in equation (2.10). If = 0, 

A|b = (2.12) 
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In time dependant flows, the quantity A is conserved along the motion; A is the inverse of the 
potential vorticity; it is also the load of Lynden-Bell and Katz (1982). The conservation of 
A is equivalent to the conservation of circulation along closed contours (Katz and Lynden- 
Bell 1985). Now consider all the closed loops of constant A that can be drawn in our 
trial configuration and parametrize them with the value of the circulation along each loop 
27rQ;(A) ; thus 

a{X) = 7^C(A) = ^ (fw-dR (2.13) 



27r ' ' 27r jx 

a{X) equals zero at the fixed point O of fig. 1 and along the boundary a{0) — aB- For 
simplicity we assumed that a increases uniformally between and as- Having defined 
a{x^) at every point of our trial flow, we now define an "angle" P by the condition that 

= Va X V/? (2.14.a) 

or 

This is a first order differential equation which defines f3 up to an arbitary function of a, 
B{a) that depends on the choice of the line P = (see figure 1). We shall define /3 = as 
the positive x axis. Following (2.13) 

a=-^(pW-dR=-^ I udxdy = ^ ! dad(3 = a{X)^ i d(3 (2.15) 
27r Jx 27r J 27r J 2tt Jx 

The domain of (3 is thus < /3 < 27r. Having now an a{x^) and a P{x^), it follows from 
(2.10) and (2. 14. a) that the velocity field is of the form 

W = aV/3 + Viy (2.16) 
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which is a Clebsch form (Lamb (1945)). It has been shown in paper I that the function v 
is single valued . The value of v is defined if we impose mass conservation , i.e. equation 
(2.2) or 

V • {aVv) = V • [a{-aVP + ffc)] (2.17) 

This is an elliptic equation for u whose solution depends on one arbitary constant only 
because (t|b = 0. We chose u{a = 0) = to make u unique. 

Thus given a mass preserving trial configuration of a stationary Maclaurin flow in 
moving coordinates, there exists a labelling a,P of the fluid elements with < a < as 
and < /5 < 2n. The labelling is isocirculational in the following sense: the circulation 
along any contour I defined by a{P) depend on a{f3) only: 



W-dR^ (ta{/3)dp (2.18) 
/ Jl 

and not on x^{a,P). Any deformed trial configuration in which fiuid elements labeled 
a, (3 with coordinates {a, P) go to points {a, P) will have the same circulation along 
the contour defined by the same function «(/?) in x^ coordinates. 

Reciprocally let x^{a, (3) define the positions of fiuid elements in a trial configuration 
labelled by a, (3 with the topology of fig.l, with X{a) given, 0<q;<Q!b, 0</3< 27r, and 
aB given as well. Then the density of the configuration a{x^) is defined by (2.11) and 
(2.14.b) as 

and a W{x^) is defined by (2.16) in which u is given by equation (2.17) and f^a = 0) = 0. 
The fiow so defined a{x^),W{x^) is both mass preserving and isocirculational. Notice 



that b and Qc appear in W through u. The three numbers 6 and Qci are, in general, related 
to the linear momentum P and the angular momentum J of the flow; P can always be 
take equal to zero: 

P = / Wad^x = (2.20.a) 
Js 

J=Tz- Rx WacPx = Jo (2.20.6) 
Js 

2.3. Fixation of Coordinates 

For deflniteness, the origin and orientation of coordinates have to be flxed. For this 
we have only the topology of the flow. We shall flx the relative position of our trial flow 
( see flg. 1 ) in the following way: The point of zero relative velocity, where a = 0, will 
be chosen as the origin of coordinates, and the x axis oriented along one of the principal 
directions of the boundary contour a = as- This flxes in general the coordinate system of 
any trial configuration . In these coordinates R{a, (3) must satisfy the following conditions: 

R{0,P) = (2.21.a) 

CJX 

(^)a=a.,/3=0 = (2.21.6) 

In addition we must have 

x{a, 0) > 0, y{a, 0) = (2.21.c) 
because we chose /3 = to be the positive x- axis. 

2.4. The Energy Criteria for Stability 

The following criteria of stability referes to what we called in paper I the strong 
energy condition, in which mass conservation always holds. Among all mass preserving 
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and isocirculational trial configurations x^{a,P), that satisfy (2.21), only those for which 
the energy 



E = 



a(fx (2.22) 



under small displacements 

R{a, (3) R{a, (3) + i{a, /5) (2.23) 

is stationary, = 0, for given P and J (AP = A J = 0), are real physical flows in that 
they satisfy Euler's equations . Physical flows, with minimum of energy are stable with 
respect to small perturbations. The condition of stability for flows with given linear and 
angular-momenta are thus 

A^E > at AE = AP = AJ = with A^P = AV ^ (2.24) 

As emphasized in paper I, (2.24) is a necessary and sufficient condition of stability if 
the gyroscopic term of the Lagrangian, the sum of those terms that are linear in time 
derivatives, is equal to zero. The explicit expression for the gyroscopic* term /S?G obtained 
in 3-D flows holds here with minor changes. Thus, the necessary and sufficient conditions 
of stability for flows with given linear and angular momentum are equation (2.24) with 

A^G = (2.25) 



3. Application to Maclaurin Disks 

* Here we deal only with the strong energy principle and drop the sub index "strong' 
in (6.9) of paper I 
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3.1. Stationary Configurations 

The following is taken from Binney and Tremaine (1987) with some changes of nota- 

— * — * 

tions. In Maclaurin disks the fluid rotates with uniform angular velocity Q, = l^Q and the 
velocity field in inertial coordinates is thus 

Wo = nx R = nR^V(p (3.1) 

an index referes to stationary quantities, (p is the polar angle, R the radial distance. 
Maclaurin disks appear the same in coordinates rotating with uniform angular velocity 
flc- In such coordinates, following equation (2.1), the relative velocity 

Uo^{n- Oe) X R (3.2a) 

and 

6o = (3.26) 

The vorticity 

a;o = Iz • rotWo = 2Q (3.3) 

The mass preserving density 



i?2 

c^o = (^c\l 1 ~ R<a, (3.4.a) 

(To = R>a. (3.4.6) 
Pressure and specific enthalpy are given by the equation of state 

P = -Ka^ (3.5.a) 

2 ^ ^ 
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h = ^na^ (3.5.6) 
The gravitational field in the plane of the disk 

$0 = l^lR', - (3.6) 

VLq is the angular velocity of a test particle on a circular orbit. Euler's equations (2.7) 
relate f2 to fio; cfc and a: 



^-" = ^1-^-^ (3.7) 



'0 2 



The mass M and the angular momentum Jq define ac and VL for any given radius a: 

2'jr 

M = —acc? (3.8.a) 

Po = (3.8.6) 
Jo = — Slcrca = — - — (3.8.C) 

iO OTT 

where 

Cb = 27ra2n (3.9) 

is the circulation along the boundary. 

3.2. (a, f3) Labelling and the Load X{a) in Maclaurin Disks. 

The labels defined in section 2 are here as follows. Lines of constant load — = const. 

(JJO 

are circles. The circulation along lines of constant load 

C{R) = 27rni?2 (3.10) 

Thus, a{x^) defined by equation (2.13) is proportional to R^: 

0<a = ^ = nR^ <aB = ^ (3.11) 
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The load X{a) is correspondingly 



A(q;) = ^ = AcW1-— (3.12.a) 
(jjQ V as 

where, with ac given in terms of M by (3. 8. a) and related to Cs in (3.9), we have 

Ac = ^ (3.12.6) 
P is defined by equation (2.14.b) which reduces here to 

If /? = is the X axis, then 

P = ^ (3.14) 

Equations (3.1), (3.11) and (3.14) give 

W = aVP (3.15.a) 



or 



Vi/Q = (3.15.6) 



and with i'o{a = 0) = 



z^o = (3.16) 



3.3. Trial Configurations of Maclaurin Flows that Differ Slightly from a Maclaurin Disk 
Consider now mass preserving and isocirculational trial configurations R{a, /3) that 
differ slightly from a Maclaurin disk Ro{a, j3) 

R{a,P) = Ro{a,P)+i{a,l3) (3.17) 
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with 

\i\ « a (3.18) 

Instead of {oi,(3), we may conveniently use {R,(f) as labels. A perturbed configuration is 
drawn in figure 1 with coordinates taken in such a way that formulas (2.21) hold. R is 
given by equation (3.17), and R as Rq satisfy both equation (2. 21. a); we must then have 
also 

f(i2 = 0, (/?) = (3.19.a) 

Equations (2.21.b and c) will be satisfied by removing from ^ any displacements that 
amounts to a change of the cut /3 = by an arbitary infinitesimal function of a i.e. 

— * 

6/3 = T{a) . In terms of equation (3.11) and (3.23) [below where 5 is related to ^ taking 
index K = 2], we may equivalently write that condition as 

R^S/S = 1^ • f X ^ = R^t{R^) (3.19.6) 

Equation (3.19.b) includes arbitrary rigid rotations, 6(3 — const. 

The density a of our Maclaurin flows change from ctq given by equation (2.19) in which 
X{a) is the load of Maclaurin disks defined by equation (3.12) with the same Xc and the 
same as- These conditions preserve indeed vortex fluxes. Thus the change in a 

Aa = a-ao = -agV ■ f (3.20) 

The velocity field W is given by (2.16) and therefore the change in W 

AW = W-Wo^ -V^- Wo + VAv (3.21) 

Ap itself is defined by the perturbed mass-preserving equation A[V ■ {crU)] — 0. To obtain 
this equation it is convenient to introduce the often used 5-variations appearing already in 
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(3.19.b), associated with changes of the functions a^{x^) that leave the hnes of constant 
fixed. Such changes are defined by 

a^{R + i} + Sa^{R + i) = a^{R) (3.22) 

To order 1, 

5a^ = -i- Va^ (3.23) 

Just hke Aa^ = 0, now 5x^ = 0. The greatest quahty of 5-variations is that they commute 
with derivatives 

(5V = VS (3.24) 

One easily find that 

5 = A - f • V (3.25) 

We shall move from A to (5 variations and back at our convenience. Thus Act in equation 
(3.20) gives 

Sa = -V- (aoi) (3.26) 

and AW in eq (3.21) gives 

5W = -Vi-W + VAu-fvW (3.27) 

Let us now go back to the equation for Au or Su. Following equation (2.2), we have 

S[V ■ {aU)]\o = V • {SalJo + aoSU) = (3.28) 

With equations (2.1) and (3. 2. a), one can see that (3.28) is an elliptic equation for Sv: 

V • {ao^Su) = -V • [SaWo + (To5{aVP)] + 5b ■ Vo-q -l^-Rx ^Saflc (3.29) 
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Equation (3.29) with Si'{a = 0) = defines Su uniquely. Su is a non homogeneous hnear 
expression in 6b and Qc- The latter will now be retrieved from the conditions that perturbed 
Maclaurin flows have same linear and angular momenta as the Maclaurin disks. 

3.4- Global Constraints 

— * 

a. Fixation of Sb 

The equation (2.2) for mass conservation implies the following identity: 

• (aU) = V • (aUR) -aU = (3.30) 
With (2.6) and (2.1), we have thus 

P = M{b + QcX Rcm) (3.31.a) 

where 

RcM = J ^^d^^ (3.31.6) 
defines the position of the center of mass. For our Maclaurin disks 

Po = bo = RcMo = (3.32) 

If we want to keep P = 0, we must take 

Sb=-ncXSRcM (3.33) 

This equation defines 5b, that appears in (3.29), in terms of I* and Qc- 

b. Fixation of Qc 

r^c is arbitary in a Maclaurin disk but as soon as we destroy axial symmetry, we can 
see the flow rotating, flc will be defined by the condition that Maclaurin flows have the 
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same angular momentum as the Maclaurin disks. Turning our attention to J, defined in 
(2.20.b), we obtain 

AJ=1,- Ji^xW + Rx AW)a(fx (3.34) 
At the stationary point (Maclaurin disks) the first term on the right hand side 

J (x Woaod^x = J ix{nx R)aod^x J i' R(^od^x (3.35) 

while, using (3.1) and (2.6), we can write the second term as 

j Rx A#|o(7od^x = y" ^ X (-Vl*- Wo)(Tq<Px = -fl J f Rao(fx (3.36) 

Since (3.35) + (3.36) = 0, we see that 

AJ|o = (3.37) 

Thus AJ|o = does not define Qc- However J|o = does. Indeed, following equation 
(3.34), 

AV|o = 1z- J{2ix AW\o + Rx A'^W\o)ao(fx = (3.38) 

which contains Av and A^v that depend both linearly on ilg- Notice, however that A'^v 
appears in the following integral 

J Rx VA'^u\oao(fx = j ^ ^ {RA'^v\Q(TQ)(fx = ^z ^ AV|o(7o^ • dR (3.39) 

The latter is zero because A'^v must be small and ctqIb = 0. Therefore Qc will only be 
defined in terms of A^J = through Au that appears in AM^|o and A^W'Iq. The explicit 
expression for Qc is somewhat complicated and will be calculated in due place. 
3.5. The Energy Principle 
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We may now apply our strongly constrained energy principle A^£' > to find stability 
limits in Maclaurin disks. With IS.v defined by (3.29) in which fig is given by (3.38) and 
bh by (3.33), we can evaluate SW given in (3.27) and (5(T in (3.26). Here is a differential 
equality convenient for calculating b?E. It is obtained from (2.22) in appendix A 

IS^E-^IS^J =8'^E-m'^J = J{a{SWf + 5a5{h+^)}\o(fx (3.40) 

4. Spherical Harmonic Decompositions of the Perturbations 

4-1- Spherical Harmonic Decompositions of ^ , Sa, Sh and 5^ 

a. ^ defined by scalar functions 

It is a good thing to define ^ in terms of two independant non dimensional infinitesimal 
scalars rj and i/j as follows: 

f = a^[V77 + rotip], -0 = l^V (4.1) 
77 is thus defined in terms of ^ by 

Ar/ = ^ V • i (4.2) 
Some boundary conditions are needed to make rj unique. If we take, say, 

v\b = (4.3.a) 

equation (4.2) has a unique solution. Equation (4.2) represents also the condition of 
integrability of (4.1), considered as a set of two first order differential equations for ip, 
given ^ and rj. Thus the i/j equations are integrable and define up to a constant. This 
constant we shall fix by asking: 

V'lc = (4.3.6) 
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b. Spherical Harmonic Decomposition of r] and ip 

We now decompose rj and in normalised spherical harmonic functions of x ^md ip 

with 



X= — = V1-^' 0<x<l. (4.4) 

We have 

oo oo 

^ = E E VimVrix)e''^'' + c.c. (4.5.a) 



l=m m=0 

oc oo 



^ = E E ^^-7'r(x)e^^'^ + C.C. (4.5.6) 

rjim and are arbitary, infinitesimal, complex numbers related by equations (4.3). As- 
sociated Legendre polynomials are defined in the range — 1 < x < 1, but rj and i/j sixe only 
defined in the range < x < 1. We may extend r], i/j to the negative region of x in any 
way we want. However, these functions have bounded gradients at x = 0. Indeed, since 

IdnTjl < oo, \dRip\ < oo (4.6. a) 

\-^d^r}\<oo, |45<pV'l<oo (4.6.6) 

Jri Jri 



hold and since 



dR^- — VT^d^, (4.7) 
ax 



then. 



d^r]\^=o = 5xV'lx=o = 0. (4.8) 
Equations (4.8) will be satisfied if we make symmetrical continuous extensions 

^(x) = V{-X), (4.9.a) 
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V'(X) = V'(-X), (4.9.6) 

whose expansions are given by (4.5) with (l — m) even. The expansions in the domain 
< X < 1, of arbitrary continuous 77 and \l) that satisfy (4.6) everywhere are then given by 
(4.5) with (l — m) even [Arfken 1985] . The boundary conditions (4.3.b) gives the following 
relations among the V'io's: 

00 

^^■0^0 = 0, / even (4. 10. a) 

while equation (4. 3. a) gives a relation among rjimS: for every m > 0: 

XI VimVn^) = 0, l-m even (4.10.6) 

l=m 

Equation (4. 10. a) defines, say, V^oo in terms of '0/7^0,0) ^ind equations (4.10.b) define, say, 
r]ii in terms of 7]i,mi^i. 

c. Spherical Harmonic Decomposition of 5a, 5h and 5^ 

We can now calculate the spherical harmonic decomposition of da with equation (3.26); 
from (3. 4. a) and (4.4) we obtain 

cro = crcX- (4.11) 
With equation (4.1) Sa can be written as: 

(5(7 = -aVc[V- (xVr/) + Vx X VV']. (4.12) 
and from equation (4.12) and (4.5) we find that 5a has the following expansion 

5a = acJ2J2'^^'^ ^'""^ + l-m even (4.13.a) 

l=m m=0 

in which 

crim = [1(1 + 1) - 'm'^]'nim + imiplm = kimVim + imipim (4.13.6) 
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With the expansion of 5a we obtain directly through (3.5.b) the spherical harmonic de- 
composition of 5h: 



5h = Znal (TlrnVrixV'^'^ + c.c. I - m even 



(4.14) 



l=m m=0 



and from the varied Poisson equation 



(4.15) 



where ^^(-z) is the Dirac function, we find the solution 5$ which has been given by Hunter 
(1963): 

oo oo 

5$ = 2$ij XI 5^ ^''"^r(x)e'""^ + C-C. I -m even (4.16.a) 



in which 



1=171 m=Q 



^Im = —gimO'lr 



(4.16.6) 



where 



(/ + m)!(/-m)! , ^ 

9lm = 22/-l[(l^),(/^)!]2 < 1 ^ > 0' ^ - ^ ^^^^ 



(4.16.C) 



— * 

4-2. Spherical Harmonic Decomposition of 5W 

Following equation (3.27) and equations (3.1) and (4.1), SW may be written as a 
gradient plus a rotational, always a convenient form for vector fields: 



5W = a'^n[VC + rot{2r]l^)] 



(4.17.a) 



in which 



C=;^(A^-e-Wo)-2^ 
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(4.17.6) 



Notice that not only is C <S 1 , |VC| must be small with respect to a ^. Au is defined by- 
equation (3.29) and therefore, with (3.29), we obtain an elliptic equation for 

V • ((ToVC) = ^77[-V(7o • rot(2r]Q) -(Q-Qc)xR- ^Sa + Vctq • Sb]. (4.18) 
2a^iz 

We now expand C in spherical harmonics 

oo oo 

C = E I] CimVr{x)e'"''' + C.C., I - m even (4.19.a) 

l=m m=0 

We also use equation (4.1) and (4.5) to find the expansion of 5b. With equation (4. 5. a) 
for 77 and equation (4.13) for Sa, we obtain, after some straightforward but a little tedious 
calculations, that 

rrrr ^ 

Clm = —r^[2mm - (1 - -^)(rim\, {l,m) 7^ (1,1) (4.19.6) 

Kim " 

and 

Cii = iWii - 2r7ii] = -i>ii - irjii- (4.19. c) 

Equations (4.19) and (4. 5. a) inserted into (4.17.a) give the spherical harmonic expansion 
of 6W. 

4-3. Spherical Harmonic decomposition of /S?E and /S?J 

A'^E is given by (3.40) with A^J — that defines fie- Near the stationary point. 



A^E 



J{a{5Wf + Sa5{h + ^)}\od^x (4.20) 



To obtain the spherical harmonic decomposition of A'^E, we replace da , 6h, 5$ and dW 
by their respective expansions (4.12) (4.14) (4.16) and (4.17), in A'^E [detailed calculations 
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are given in appendix B]. The result is as the follows: 

°® 2 2 o 02 

/^^E = A^ac^'^a' J] A{hm - ^)\mn.\' + [^(1 - - 1 + (1 - ^.^)^]k/^| 



, „ ""I'm "'Irn ^' ^' 

(<=m)>l,m=0 



(4.21) 



The parenthesis (1 — ^) is given by J = (the calculation is in appendix C): 



^-'ir = —^ ^ I '"'^12 — (m^O),(/,m^l,l) (4.22) 



Fictitious Marginal Instabilities 

Direct inspection of A^E in equation (4.21) reveals that the following modes are 
absent: rju, an = rjn + iipn and ipio- Equation (4. 10. a) defines ■0OO) equations (4.10.b) 
define rju. The other modes, •0/o I > and an, must thus represent symmetries. It is 
important to figure out what these are. 

The absence of an = ijn + iijjn in A'^E is related to A^^'s invariance for uniform 
translations of the coordinate axis. Translational invariance is disposed ofi^ by imposing 
equation (3. 19. a) whose spherical harmonic decomposition is 

J2Kl + '^)ivii + iM = Q (4.23) 

From which we see that an — rjn + ii^ii can be retrieved in terms of rjn + iipn for I ^ 1, 
Hence, the origin of coordinates is redefined, in accordance with condition (2. 21. a). 

The coefficients il^io, I > are defined by equation (3.19.b). If we use (4.1) and the 
spherical harmonic decomposition (4.5), we obtain from (3.19.b) 
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Terms dependent on (/? have been removed from the left because r is independant of ^. 
Equation (4.24) depend on t/j/o, I > only; the latter can be retrieved explicitly with two 
obvious integrations: 

V^/o = /' ri{x)[ f r(x')x'c?xVx I > (4.25) 
Jo Jx 

The way to keep the cut /3 = fixed on the x axis is simply to set tpio = for i > 0. Notice 
that '020 = alone avoids rigid rotations. 

5. Necessary and Sufficient Conditions for Stability of Maclaurin Disks 

5.1. The Gyroscopic Term 

Following section 2.4 the inequality A'^E > 0, calculable from (4.21) with (4.22), gives 
sufficient conditions of stability. A'^E has the form ^^A^ — B > 0. Ifi?<0, the disk 
is stable for any For instance, Maclaurin disks are stable to any perturbation that 
keeps the same densities at the displaced points {6a — 0). We are naturally interested in 
perturbations that might upset stability and for which B > 0. If B > then must 
satisfy the following inequality 

Q—^^ ^{l=m)>l,m=o(^ ~ 9lm)Wl7n\'^ 

(5.1) 

The condition (5.1) becomes suflftcient and necessary when the gyroscopic term for dynam- 
ical perturbations A^G is zero (see paper I). It is therefore important to calculate A'^G. 
For dynamical perturbations, is a function of the time t as well as of ct, /3 . Gyroscopic 
terms are those bilinear functionals of ^ and ^ = §f which appear in the Lagrangian of 
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the dynamically perturbed equations. We have calculated /S?G in paper I. For perturbed 
Maclaurin Disks, A^G of paper I reduces here to a rather simple expression [see appendix 

D] 

A^G = 2 ^ j i- AWao(fxdt (5.2) 

We can make a spherical harmonic decomposition of A^G, with ^ given by (4.1) and (4.5) 
and AW given in appendix by (C.5) . The only novelty in these calculations is that rjim 

— * — * 

and ipim are now functions of t. A straightforward subtitution of ^ and AW in equation 
(5.2) leads to the following expression for A^G: 

A^G = AiTrac f dt V mO{( ^^~^^ - -^)at^aira+ 

Jtn 1 ^1 1 '^Im 



*o /,m7^1,l 

kim 



(1 ~ -^)Vim'nim - 2(t^^ -^)'nim'^im} + couiplex Conjugate (5.3) 



in which (1 — ^) has to be replaced by the value given in equation (4.22). Necessary and 
sufficient conditions of stability are obtained from (5.1) when A^G — 0. 

5.2. Symmetric and Antisymmetric Single-Mode Perturbations 

For symmetrical and antisymmetrical modes, we take either real or imaginary com- 
ponents of Tjim and V^/rn in the spherical harmonic expansion, so the Fourier expansion 
contains either cos{mifi) or sinirrup). For such perturbations, the written part of A^G 
is imaginary; adding the complex conjugate makes thus A^G = 0. For one single mode 
{l,m), equation (5.1) reduces to: 

Q<Qlm = , ^5 4) 
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in which 

1-^ = '-^ (m^O),(/,m^l,l) (5.5) 

kim 

It is advatageous to introduce a single arbitary variable z and a constant x to analyze the 
lower bound of Qim- 

z = —kim-^^ — oo < ^ < oo (5. 6. a) 



and 



x= 0<x<l. (5.6.6) 



Im 



In term of z 

/n (-'- ~ 9lm) /r T ^ 

= -Fk(^ ^'■^■''^ 
where ^4(2;) is a polynomial of order 4 in z: 

P4{z) = 1 - 7:^{(- - 1)[(1 - x)z'^ + 4(1 - x)z^ + + 4z] + -} (5.7.6) 

P4(z) has one and only one real maximum for any x (see figure 2). We are only interested 

in values of z for which P4{z) > 0. For P4{z) < 0, A'^E > for any value of Q. The 
maximum of P^iz) is obtained for: 



Zmax{x) = -^^[{VT^ - 1)^ + {VT^ + 1)^] - 1 (5.8.a) 



for which 



i"4(^)maa. = l-T:^{(^-l)[(l-a;)4a.+4(l-a:)2;i,, + 64,,+4w] + ^} = 1-^ 

(5.8.6) 

26 



The function y{x) appears in figure 3. To any pair of values {l,m) corresponds a value x 
defined by (5.6.b) and a point {x,y) on the curve. Following (5. 8. a) and (5.4), we must 
thus have 

Q ^ {Qlm)Tnin — i/Or)~ i^'^) 

The smallest minimum of Qim is obtained for {l,m) = (2,2) for which (Q22)mm = |- 
Therefore the necessary and sufficient condition of stability with respect to symmetric or 
antisymmetric single mode perturbations of Maclaurin disks is 

Q<1 (5.10) 

5.3 Comments 

a) Binney and Tremaine have given the following dispersion relation for the cu - modes of 
dynamical perturbations: 

UJ^-UJr{^n'' + kirn[^l{l-girn)-^^]} + 2mn[nl{l-girn)-n^] =0, Ur = U-TuVL (5.11) 

Stability holds if the non spurious o^^'s are real roots. Equation (5.11) is a third order 
polynome. The condition for a polynome of the form x^ + a2X + 03 = to have only real 
roots is given in standard handbooks [for instance: Schaum's Mathematical Handbook 
(1968)] 

(f)^+(f)^<0 (5.12) 
Inequality (5.12) for equation (5.11) is exactly our inequality (5.9). 

b) For I — m modes , A^G — for asymmetrical modes, that is, for an complex [ iju does 
not appear in A^G nor in A'^E] . However, the stability limit is still given by equation 
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(5.9) with I = m, for which y = 1 and kim = I' 



Q<Qii = yrry' = ^21^ (^-i^) 

c) For radial modes, m = 0, A'^E does not depend on fie, A^J = and A^G = 0. 
Moreover, since [see equation (4.13.b)] 

aio = 1(1 + l)mo (5.14) 

all Z-modes are decoupled, A'^E is a sum of squares of independent modes. In this case 
one obtains stability limits for non single modes. The most unfavorable limit, however, is 
again that given by Qio for which j/ = 4 and ku = 1(1 + 1): 

Q<Qio = r^^^ a>2), 910 = 22/-r|(|),]4 (^.15) 

d) For coupled I = m modes, A'^G and we obtain only sufficient conditions of stability. 
In particular if we couple a pair of modes I = m and I' = m' , sufficient conditions of 
stability go from Q < Qu for Z' = Z to 

Q<l-gi'i' (5.16) 

for I I'. The smallest values of (5.16) is | which is also the secular limit of stability 
found by Binney and Tremaine (1987). 

Joseph Katz's contribution to this work started during his stay at the National As- 
tronomical Observatory in Mitaka, whose hospitality he gratefully acknowledged. 
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Appendix A: lS?E - QA'^J Variational Identity (3.40) 

We start from (2.22). Varying (2.22), using (2.8) for e and (2.9) for $, and calculating at 
AE = 0, we have 

A^E^ J{{AWf + W ■ A^W + A[C-V{h + ^)]}\o(7od^x {A.l) 

By substracting from A'^E, QA^J where A^J is given in (3.38), we obtain 

A'^E-nA'^J = J {(A#)2 + #-A2#-li-(2fxA# + ^xA2#)+A[fv(/i+$)]}|o(7oc^^2; 

(^.2) 

The A'^W drops out of equation (A. 2) because of (3.1) . The terms containing AW can be 
written as a difference of squares: 

{AWf-rii2ixAW) = {AW-flxif-{Clxif = {AW-i-VWof-fl^i^ = {swf-n^^ 

{A.3) 

In the last equality we have used (3.25). So equation (A. 2) with (A. 3) can now be written: 

A^E-nA'^J = J{{6Wf -n^i^ + A[i-V{h + ^)]}\oao(fx (AA) 
With the following operator identity 

AV = VA - • V (A.5) 
and with equation (3.25), we may rewrite the /i + $ term in (A.l) as follows: 

A[f-V(/i + $)] =f-VA(/i + $)-f-V^V(/i + $) =^W(/i + $) + fvV(/i + $)-|' (A.6) 
But from Euler's equations (2.7) we see that 

V(/i + $)|o = -#o • VWo {A.7) 
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Subtituting (A. 7) into (A. 6) and A[( ■ V(/i + $)] back in (A.4), we find that il'^^ cancels out. 
Finally with Sa given in (3.26), and with some integration by parts, we obtain the following 
much simpler form for A'^E — flA^J which is that written in (3.40): 

A^E-nA^J^ J{a{SWf + 6aS{h + ^)}\od^x (A.8) 

Appendix B: Spherical Harmonic Decomposition of A'^E 

Let us start from (A. 8) and take J = 0. We shall set 

A'^E = A'^Ek + A'^Ep (B.l.a) 

in which 

A^Ek = J[a{5W)Vx (B.l.b) 

and 

A'^Ep = J[da{dh + d^)]\o(fx. {B.l.c) 

Inserting (4. 17. a) in (B.l.b) gives: 

A'^Ek = a^n^ J {{VCf + 2V • [Crot{2r]l^)] + [rot{2rjl,)f}ao(fx {B.2) 

Using identity (2.6), we can integrate (B.2) by part, and with (To|b = 0, obtain: 

A'^Ek = a'^n'^ J [-CV • (aoVC) - 2CV(To • rot{2r]l,) - 4r^V ■ iaoVr])](fx. {B.3) 

The spherical harmonic decomposition of A'^Ek is obtained by replacing in equation (B.3), C 
by (4.19), ctq by acX (equation (4.11)) and r] by (4. 5. a). The result comes out as follows: 

oo 2 2 o 

A'E,=47racn'a' J2 [4(fcz.n - ^)|r/.^P + ^(1 - ^)'k^n^P] (BA) 

l=m,m=Q 
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Notice that these are aheady some of the terms of /S?E given in (4.21) 

The spherical harmonic decomposition of A'^Ep follows by inserting in (B.l.c) the respec- 
tive expansions of 5a in (4.13), 5h in (4.14) and 5^ in (4.15). For A'^Ep, we obtain : 



oo oo 



A^Ep = 2na^Ql J] ^ 'Jimi^L + Skctc^L] + c.c. (5.5) 

l=m Tn=0 

With (4.16.b) and (3.7), A'^Ep can also be written 



A^Ep = 47raca''n^ J2 h'^ + - 9im)^]Wim\'' {B.Q) 

l=m,m=0 

The sum of A'^E^ given by (B.4) and A'^Ep of (B.6) is the expression of A'^E written in (4.21). 

Appendix C: Spherical Harmonic Decomposition of A^J 

Following (3.38): 

A2j= 1^- J {2ix AW + R X A'^W)\oao(fx. (C.l) 
By varying AW in (3.21), we obtain A'^W: 

A2#= -2V|'- A# + VAV (C.2) 

Inserting this into (C.l) gives: 

AV = 1^ • y"(2f X AW -2Rx • AW)\oaod^x. {C.3) 

Notice that VA'^u do not contribute to (C.3) [see equation (3.39)]. Following equation (3.25) 
we can write: 

AW = SW + i- VWo (C.4) 
Using (3.1) for Wq, (4.1) for I*, and (4.17.a) for 5W, we obtain for (C.4) 

A;^ = a^il[ViC + V) + rot(r]l^)] (C.5) 
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Inserting (C.5) and (4.1) into (C.3) gives 

J = 2a^m^ • y"{[V?7 + rotif] x [V(C + V') + rot{r]l^)] 
-Rx V[Vr7 + rotip] ■ [V(C + V') + rot{r]l:,)]}ao(fx. (C.6) 

The spherical harmonic decomposition of A^J is obtained by replacing in equation (C.6) ( by 
(4.19), ctq by acX (equation (4.11)) and rj, ip by (4.5). Using also (4.13.b) for aim, the result 
comes out as follows: 



AV = Anaca'n J^^iJ^rn - m')\vim - + (t^ - (1 - + cc] 

l^lm Kilm i ^ Kim 

C^Zm |2 , I |2 /I Oc^Y^m^ 



f^im l^im '^Im 

and demanding A^ J = we see that: 



Im 



^-'-7^ = ^^ ^^W^ni — (m^O),(/,m^l,l) (C.S) 

which is equation (4.22). 

Appendix D: A^G For Maclaurin Disks 

Equation (6.9) of paper I reads as follows (we drop the index strong): 




A^G = j J [VAi/D ■ VAi/g - ^ ■ • Wo]aod^xdt 

[AffcD ■ VAus + Affcs • VAuD + AffcD -Vi-Wo 
+ (AQcD xWq + QcX VAi/£)) ■ i]\oao(fxdt, (D.l) 

u and Au are defined by the time dependent equation of mass conservation and it's first 
variation. Av may be decomposed into a "dynamical contribution" Aud which is zero when 
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^ = and a "steady" part Avs which is the one we encountered in perturbed steady flows; 
Avs is the Au = 5v of this paper which appears for the first time in (3.21); here, however: 

{Av)paperl = Avd + Avg {D.2) 

The equation of mass conservation is what defines equations for Au^ and Aus'. 

V • {aoVAvD) = V • [cTo(f + A?7cc)] {D2.a) 

V • {ao^Aus) = V • [ao{Uo ■ Vf + • Uo) + ^oVf • rfe + Aff.s)] (D.S.b) 
Similarly, ffc and At/c (equation (2.1.b)) have a dynamical and steady contribution and we write 

Affc = AffcD + Affcs (DA.a) 

AffcD = Abo + AdcD X R (DA.b) 
Affcs = Abs + ^cs X i+ Ajlcs X R (DA.c) 
The four last terms of (D.l) can be written as follows with the help of 

APd = J VAuDcro(fx (D.b.a) 

and 

AJd = Tz- J VAuDaod^x (D.b.b) 

and APs and AJg which are the same as our AP leading to (3.33) and A J given in (3.34); 
now using the definitions in (D.4) we see that the four last terms of A'^G (see D.l) 

- J J [Affcs ■ VAuD + AffcD ■ Vi- Wo 

+ (AilcD X Wo + ^^c X VAi/fl) • ^ \oao(fxdt, (D.6.a) 
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become: 




[ J AffcD ■ VAus(rod^x]+ 
Abs ■ APd + AQcs ■ AJd - Abo ■ APs - AQcd ■ AJsdt, {D.G.b) 
Inserting (D.6.b) back into (D.l) gives thus: 

A^G = 2 y y [^Ai^D ■ VAi/5 - i- • Wo]aod^xdt - jiij '^^VcD ■ ^Aiys(Tod'^x] + 
Abs ■ ^Pd + Aftcs ■ AJd - Abo ■ APg - AQ^d ■ AJsjdt, {D.7) 
Conservation of linear and angular momentum works now as follows. AJs and AJd vanish 

— * — * 

identically due to the symmetry of the problem. Define Abs and Abo by setting separately 
APs = and APd = 0. With this, A^G reduces to 



A^G = 2 y y [VAuD ■ VAus -^-V^-Wq- AffcD ■ V5u]aod^xdt {D.8) 

Integrating the first and the last term of (D.8) by parts: and replacing finally Aus by our 
present Au = Sv defined in (3.29), we obtain: 

A^G = 2 y y {Ai^V • [doiAffcD - VAi^d)] - f - V^ - Wo(To}d'^xdt {D.9) 

and with (D.3.a), (D.9) can also be written: 

A^G = 2 y y [-AzvV • (aoi) - f - vf - Wo(To]d^xdt (D.IO) 

And again, integrating by parts the first term of the integrant we obtain: 

A^G = 2 y J i- {VAu - VC • Wo)aod'^xdt (D.ll) 

and finally using (3.21) in (D.12), we get the simple form of A^G written in equation (5.2): 

A^G = 2 [ [ i- AWaod^xdt {D.12) 

Jtn J 
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Figure Captions 

Figure 1: Examples of a trial configuration. The drawing shows some a — const and P = const 
lines and the positioning of coordinate axis. 

Figure 2: This represents one P^iz) > ioi x = ^ coresponding to {l,m) = (4,2) and a 
Zmax — —0.734. The hatched part of this curve corresponds to perturbations stable for any Q. 
Figure 3: (No Caption). 



36 



